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PREFACE 


This is the final report on "Supplemental Selected 
Studies of VHF/UHF Communications for Planetary (Mars, 
Venus) Relay Links", prepared by the Astro-Electronics 
Division of the Radio Corporation of America for the Ames 
Research Center, National Aeronautics and Space Adminis - 
tration under Contract No. NAS2-3772. The report covers 
work accomplished during the period from March 17 to 
September 30, 1967 . All of the information presented here 
was initially highlighted in monthly reports and discussed in 
detail at the final oral report presented at the Ames Research 
Center in August 19617. Copies of the presentation material 
used in the oral report have been delivered to the Ames 
Research Center. 

The major contributors to this report are Dr. H. Staras, 

Mr. S. H. Koth, Mr. W. Bisignani, and Mr. H. G. Rouland. 

Dr. Staras performed the multipath analysis for an exponen- 
tial-type surface model (Task IA) , and Mr. S. H, Roth and 
Mr. W. Bisignani were responsible for the performance of 
the fast-fading experimental investigation of a wideband non- - 
coherent FSK system (Tasks 3, 4, and 6). Mr. Rouland provided 
the program technical management. 
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SECTION I 


INTRODUCTION AND SUMMARY 


A. INTRODUCTION 
1. Background 

The specific studies performed by RCA during this program supplemented 
and extended available analytical information relative to communications problems 
between a probe entering a planetary atmosphere and a nearby spacecraft bus, 
based on two probe trajectories representative of entry from a flyby mission and 
from an orbit mission. The study comprised four tasks: 

© Investigation of multipath environment using an exponential surface 
model of Mars and Venus; 

• Experimental determination of the effects of fast fading on the wideband 
noncoherent FSK relay link with hardwire bit synchronization; 

© Experimental determination of the effects of fast fading on the relay- 
link receiver with bit synchronization; and 

© Experimental determination of fast fading on the distribution of errors 
as a function of the mission profile. 

The original multipath analysis performed under this contract was based on an 
analytic surface model (two-dimensional Gaussian distribution) and led to useful 
results in defining the multipath environment (average scattered power and fading 
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bandwidth from an extended surface)- Radar data from the moon and other planets 
indicates that the planet surface may be slightly rougher or more diffuse than the 
analytic surface model implies. At the beginning of the study, it was felt that a 
non-analytic surface with an exponential correlation function would be a more 
representative surface model than the analytic surface model. 

The purpose of Task 1A, "Multipath Propagation Study, " was to refine the multi- 
path analysis by utilizing the non-analytic surface. The scattering function for 
an extended surface could not be solved in closed form but was solved using a 
numerical procedure. The scattering formula derived for an extended rough surface 
with an exponential correlation function yielded results of such doubtful accuracy 
as to make the validity of the approach highly suspect. Analysis of the basic 
formulation showed the conservation of incident power to aggregate of the re- 
flected power is not satisfied for differential areas with unity reflection co- 
efficients. Attempts to "doctor" the formulation were not successful. 

The basic analysis of the scattered power from an extended exponential surface 
is presented in Section If, "Multipath Propagation Study (Task 1A). "Additional 
discussion of this analysis, including recommendations, is presented in Section 
I, Paragraph B, "Special Notes on Multipath Propagation Study (Task 1A)„" 

The previous study considered the effects of slow fading, which is the worst case, 
on the performance of & wirieL<>'- d noncoherent FSK lini . For the low data 
rates of interest (100 to 300 bits per second), the multipath analysis showed that 
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the link would be operating in a fast-fading environment. In order to experimentally 


determine the effects of fast fading, a simulator was assembled with three discrete 

n 

fading bandwidths, continuously variable y * and doppler offset of the reflected 
signal with respect to the direct signal. The simulator discription and the effects 
of fast fading on the wideband noncoherent FSK system are discussed in Task 3, 
"Experimental Investigation of Fast Fading. " The effects of fast fading on bit 
synchronization is discussed in Task 4. 

In the case of slow fading, the fading period is long compared to a data bit. When 
the link is in a fade, many bits are affected and the errors are not independent. 
Therefore, the mean probability of error may not be as meaningful as the 
percentile probability of error, which was discussed in the previous study. The 
purpose of Task 6, "Fink Performance Criteria" is to experimentally determine 
the effects of moderately fast fading on the distribution of errors as a function of 
fading environment corresponding to three different points in the probe after entry 
into the atmosphere (near end of blackout, midway to impact, and near terminal 
impact). This is especially important since the aggregate of all the data trans- 
mitted from the probe may not be sufficient for a good statistical sample. 

2. Summary of Task Description and Results 

The study comprised analysis of multipath propagation for a non-analytic 
planetary surface model and experimental investigation to determine the effects of 
fast fading, which is expectea lo ’ a probe entry mission, on the communication 
link. A summary of results for the multipath progagation task is presented in 
Paragraph B of this section, 
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The experiments utilized the breadboard system developed earlier under this con- 
tract. Unless otherwise specified, the receiver was operated at a data rate of 
lOObits per second which corresponds to a WT of 100. As part of the experimental 
effort, a fast fade generator was developed with the following characteristics: 

Fading Bandwidth 200 Hz, 1000 Hz and 3000 Hz 

y 2 (variable) 0 to °o 

Doppler Offset 0 Hz to 10,000 Hz 

For the above fading bandwidths and a data rate of lOObits per second, the fading 
rates (BT) are 2, 10, and 30. 

The results of the experimental determination of fast-fading environment on a 
wideband noncoherent FSK system are: 

e The effect of increasing fading bandwidth (for a constant reflected 

power) is to improve the link performance. The performance can be 
superior to the no-multipath case, but improvement beyond the no- 
multipath case cannot be utilized in the mission since it relies on 
the presence of the reflected signal, which can only be postulated, 
but not guaranteed. 

® The effect of increasing differential doppler offset (for a constant 
reflected power) is to improve link performance. 

e The bit synchronizer introduces approximately 0.5 db to 1.0 db 

degradatioi in link performance in a nonfading environment. While 
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previous tests indicated that slow fading adversely affected the per- 
formance of the bit synchronizer, the results here show that fast fading 
does not measurably degrade the performance of the bit synchronizer. 

e Fast fading (BT S 2 ) does not cause any measurable clustering of 
errors. Since the errors are essentially independent, the mean or 
average probability of error is a valid measure of link performance. 

® The limit of improvement in link performance, due to fast fading, 

increasing doppler offset, or any combination of the two, is to consider 
the received signal energy as the sum of the direct and reflected energy. 


B. SPECIAL NOTES ON THE MULTIPATH PROPAGATION STUDY 
1. Background 

In the previous study performed under this contract,' formulas were derived 
for the amplitude of the reflected signal relative to that of the transmitted (direct) 
signal, the polarization of the reflected signal, and the fading bandwidth of the 
reflected signal. The analysis was based oil an analytic planetary surface model 
which was developed initially by Beckmann and Spizzichino' for differential surface 
areas. The original effort by RCA was to determine the amplitude and distribution 
of the reflected signal from an extended surface (based on the Beckmann and 
Spizzichino formulations), which would be illuminated by a probe or satellite near 
the planet's surface. The principal assumptions associated with this approsch are 
summarized below: 
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( 1 ) The rough surface is Gaussian, with <j/X » 1, where a is the 
rms "bump" height, and X is the radio wavelength. 

(2) 'The rms slope of the surface is fairly small. Mathematically, 
this is written as a/ T « 1, where T is the correlation-length of 
the surface undulations. 

(3) The rough surface correlation function is analytic. The mathe- 
matical implications of this are discussed by Beckmann 3 , as 
well as others, 

( 4 ) Shadowing of one part of the surface by another, as well as mul- 
tiple scattering are not important considerations. This is a satis- 
factory approximation only as long as (2) and (3), above, are 
satisfied. 

The generally accepted solution for the scattered signal from the analytic surface 
based on these assumptions is that the energy is concentrated in a narrow beam 
centered about the specular direction. Review of radar data from the moon, how- 
ever, suggests that some of the energy should appear in sidelobes (implying a 
slightly rougher or more diffuse surface than is indicated by the analytic surface 
model), and that a non-analytic surface (specifically one with an exponential 
correlation function) would be more representative, More recently, another 
paper 4 appeared describing model studies of acoustic scattering from a rough 
surface which also suggested that the exponential correlation function provides a 
more reasonable fit to the experimental data than the analytic function. 
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The purpose of the study reported on here was to determine the significant multi- 


path parameters from an extended exponential surface to obtain the second order 
(sidelobe) effects associated with a surface which scatters somewhat more diffusely 
than one characterized by an analytic surface. 

2. Summary of Results 

For the analytic surface previously studied, the rough planetary surface was 
characterized by T/a, which represents the reciprocal of the rms slope of the 
surface undulations. For the non-analytic surface postulated for this portion of 
the study, the rough surface was characterized by a parameter, p, which is pro- 
portional to ( T/a) (A/a)= While T/a is large, X/a is small for a rough surface. 
The surface of Mars, Venus, and the moon are probably characterized by p's of 
the order of 1. 

For the p*s of interest, the scattering formula derived for an extended rough sur- 
face with an exponential correlation function yielded results of such doubtful 
accuracy as to make the validity of the approach highly suspect. Analysis of the 
basic formulation showed that for small values of p, the required conservation of 
incident power in the aggregate of the reflected power is not satisfied even for 
differential areas and unity reflection coefficients,. Attempts to "doctor" the 
formulation were not successful, Only for values of p greater than 10 (which are 
not of interest in this study) did the scattering form-ula appear to approach validity. 
Inspection of the numerical computations indicate that the formulations attribute 
excessive scattering power near grazing angles. 

I-? 


> 



The exponential correlation function was known at the outset of the study to be un- 
realistic 5 , since the first derivative of the function is not zero at the origin (a neces- 
sary condition for complete validity). However, it was hoped that it would furnish 
information regarding reflection distributions away from the origin ( sidelobes). In 
particular, the exponential correlation function implies that the slopes of the surface 
undulations can be extremely high, and such high slopes explain why the resulting 
scattering function is highly diffuse. However, the existence of such slopes conflict 
with a basic assumption used in deriving the scattering function: namely, that the 
rms slope is small. Furthermore, the steep slopes imply that significant shadowing 
and multiple scattering occur, which also conflicts with an explicit assumption used 
in deriving the scattering function (see Paragraph 1, above). 

It is of interest to note that the analytic surface model studied previously does 
satisfy conservation of power requirements. It is felt that a modification of the 
generally accepted solution for the reflected signal from an analytic surface will 
yield better estimates for the scattering pattern away from the specular direction. 

3. Recommendations 

The two-part analytic study performed under this contract has highlighted 
a topic in this subject of radio-wave scattering from a rough surface which deserves 
further elucidation. 

The analytic correlation function assumed in the first part of the study gave rise 
to a scattering function which was highly concentrated around the specular direction 
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and decayed very rapidly away from the main lobe. This decay was more rapid 
than experimental data seem to suggest. To accommodate the experimentally 
observed slower decay, various investigators have suggested the use of an expo- 
nential correlation function (which is, of course, not analytic at the origin). 

While the exponential correlation function does indeed give rise to a scattering 
function which decays more slowly away from the main lobe, the present study 
has demonstrated that it has another serious flaw; namely, it does not satisfy 
the conservation of power requirement. Specifically, what is meant by this is; 
the scattering function that results from the use of the exponential correlation is 
of such a mathematical form that the total power (i.e. , over all directions) scat- 
tered from an element of surface is run, equal to the incident power intercepted 
by the element of surface. This failure to satisfy the conservation of power re- 
quirement implies that its use in calculating the average scattered power from a~ 
extended surface is of doubtful validity. It, therefore, appears desirable to find 
a more satisfactory analytic formulation of the problem. The following mathe- 
matical considerations suggest that this may be possible. 

Beckmann and Spizzichino have shown that the scattering function, P , can be 
written 

P » F 2 f t J (v r) e“ g ^~ C(T) ] dr (1-1) 

s'Joxy 

where 

l+cosf^cosfl - sin 8 sin cosy 

F = 

cos0 +cos0 2 

v 2 = k 2 (sir?Q + sdi?Q n - 2 sine sine cosy) 
xy 12 1 2 
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2 2 2 
g = k a (cos0 1 + cos0 2 ) 

C(t) = surface correlation function 

Since g is assumed to be very large, it has always been assumed that C(T) need 

be known only in the neighborhood of T = 0 where C ~ 1. Therefore, one of two 

2 

T T 

assumptions have been made in the past. Either C(T) ® 1 - or C (T).« — , 

T 2 T 

the firstbeing the approximate form for an analytic correlation function, while 
the second is the appropriate approximation for the exponential correlation function. 
Upon further reflection, it appears that the analytic approximation mentioned above 
is probably valid when g » v T (namely, not too far from the specular direc- 
tion); however, it is probably not valid when g « v T (namely, quite far from 

xy 

the specular direction). It would appear that further study of the integral in 
Equation I— 1 , especially from the viewpoint of finding an asymptotic form for 
P when g « v T even if C (T) is analytic, would help clarify the whole question 
of scattering from extended surfaces. 
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SECTION n 


MULTIPATH PROPAGATION STUDY (TASK 1A) 

A SUMMARY 

This study considers the effects of a rough planetary surface on certain 
communications system parameters of interest that arise in various planetary 
exploration missions. The study represents an extension of an earlier study 
performed under this contract. It differs from the earlier one in that the cor- 
relation function that describes the surface roughness was assumed in the pres- 
ent study to be exponential rather than analytic," The physical implication of 
this new assumption about the correlation function is that the pattern of inco- 
herent scatter from the rough surface is more diffuse than it was under the as- 
sumption of the earlier study. 

It was hoped to evaluate analytically the average scattered power and the fading 
rate. Unfortunately, this hope was not realized. The diffuse scattering function 
is so complicated that it was impossible to find a satisfactory analytic evaluation 
of the total scattered power, let alone the fading rate. It was decided, therefore, 
to put the integral for the average scattered power on a computer. This numeri- 
cal computation procedure brought to light a basic defect in the formulation of 


The exponential correlation function is, of course, not analytic at the origin, 
since it is assumed f here is an exponential decay on both sides of the origin. 
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the problem, a defect that makes evaluation of the average scattered power from 
an extended surface somewhat suspect in the case of diffuse scattering, This 
vnU be discussed in some detail in the body of this section. 


Since no analytical procedure was found that would permit closed-form evalua- 
tion of the average scatter ad power, it was not possible to carry out the other 
integrations that would have yielded the fading rate. However, the basic physical 
principle suggests that the fading rate in the case of very diffuse scattering would 
be higher than in the case of directional scattering, since the size of the effective 
scattering region is larger. 


B. THE AVERAGE SCATTERED POWER 
1. Analysis 

a. Basic Formulation 

In the previous study' (i.e. , based on the assumption of an analytic 
correlation function), it was shown that the ratio of the average scattered power 
to the direct power can be written 


< P. > 
s 





(n-i) 


p 

o 



2 

— r2 - eo t: ¥po 
— 2 2 

r 2 cos (3 


exp 


tan^g 

tan^g 

o 


(H-2) 



where tan ^ ("'-' rms slope of surface undulations), 

2 2 

2 sin Oj + sin 0 2 - 2sin0j sin 02 coscp . 
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tan f3 = 


(cos 0^ + cos d^y 


and R is the Fresnel reflection coefficient for a plan surface. 

The various angles and distances defined above are the same as those defined in 
earlier study but are presented in Figure 11— 1 for convenience. In the present 
study (i.e. , based on the assumption of an exponential correlation function), Equa- 
tion (II— 1) still applies, but Equation (II-2) must be rewritten to read 


P - 
o 


^|2 ^2 J\+cos0^cos02 " sine^sin02coscpj'‘ 


,,2 4 2 2 

2irk 0 x i v 2 6 

(cos0^+cos0 2 ) 


2 

T 2 

14 = 

,2 4 r 
. k a 


3/2 


(II— 3) 


where 


2 2 

2 sin ®-j. + s *- n ®2 ~ ^ S ^ n s ^ n C ° S 




(cos 0^ + COS e,) 


4 


k = 

A C 


Equation (11-3) differs from Equation (11-2) in two important respects. First of 
all. Equation (II— 2) is frequency independent, while Equation (II— 3) is frequence 
dependent. Secondly, if T/cr»l, Equation (II-2) decays more rapidly away from 
the specular point (where tan /3 = 0) than does Equation (11-3). There are two 
reasons for the less rapid decay cf Equation (11-3) compared to Equation (II-2). 
One is that the important, parameter is not T/cr (which is very large) but rather 


T/cr -*■ kcr, where kcr is also very large. It is possible that is not much larger 




Figure n-1. Geometrical Parameters Associated with 
Scattering from a Rough Surface 


than unity, and is perhaps substantially smaller. The other reason is the presence 
cf an inverse 3/2 power decay from the specular point, which is of course a far 
slower decay than the exponential function in' Equation (II— 2) 


After many attempts to find approximate analytic solutions for Equation (ITT) using 
Equation (11-3) as the integral, it was decided that none of these were sufficiently 
reliable. Not being able to evaluate Equation (Il-I)analytically, it was decided 
to use numerical methods. Luckily, the integration over cp can be done in closed 
form to obtain a complicated sum of complete elliptic integrals. Minor manipu- 
lation of Equation (II— 3 ) permits the representation 

tj 2 |jL+cos(0 ^“0 2 " ,2s in(?i sine ^cos*" ^ 


;° * . , 
■T/H 

o o r 




4 2 2 2. 

( cos 0 ^+cos 02 ) p2(sine _j+sin© 2 ' 2 " 4p sine ^sine 2 Cos g 


(II-5) 

pdpdq 

3/2 


Substituting \|r = — _ ^ yields 


I = 


2jl. 


pdp 


o r 


4 2 2 1 

(cos0^+cos02) + P (sinG^+sin©^) ■ 


,1 3/ 2 


jt/2 

/ 


Mil 

D(^) 


di]r (II-6) 


where 


N 


D = 


l+cos(0^-02) 

L J 


- 4 


I+cos(0^-02) 


. • 2 2 2 4 

sin 02 Sin 02 sinAjf+ 4 sin 0 ^sin 62 s ^ n / 


i - 


4p2sin0^ sin 02 
k 4.2 


(cos0 -j+cc^A^) + p (sin 0 ^+sin 02 ) 


9 s in 2 i}r 


3/2 
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The integral over \k is now in the standard form of the elliptic integrals. Intro- 
ducing q 2 = (cosfy+eos^) 4 + P 4 (sinf^+sin^) 2 and 


9 4d 2 — 

k 2 = — ^ sin0 1 sin0 2 , we obtain 
ft 


jt/2 


tc/2 


n 


f |r| 2 p dp j f f+cos( 0 1 - 0 2 )] 

J , 2 J I * 2 ' 2 3/2 


r q 


o (1-k sin t}t) 


cty - 4sin0^sin0 2 |l+cos(0.^-0 2 ) 


[ 2 

I - sin — t j 

J n n 


dx 


o • (l-k 2 sin 2 ijr) 


k/2 


2 2 

+ 4sin e^sin 0^ 


/ 


• 4. ,, 

sm \lr dijr 


(l-k 2 sin 2 \j/) 


3/2 


aid 


r j 

\r\ 

|2 A 

1 p dp 

| 

/ j 

> 3 < 



[l4cos(0 1 -0 2 )] 


-1-k 


E(k 2 ) - 4sin0 sin0 [l+cosC©., -0, )1 

L 2 J k z (i-k 2 ) 


+ 4sin 2 0 lS in 2 0 9 ^ 2 ) , E„<k 2 )-2(l-k 2 ) K(k 2 ) 
1 2 (1-k 2 ) k 4 


(n-7) 


Equation (II— 7) was evaluated numerically with the aid of the following trans- 
formation (refer to Figure 11-2): 


x = COS0^ 


h 

r 


dx = -sin 0 ^d 0 ^ 


p = htan0^ 


dp -- k 2~ de i 
cos e ^ 


pdp dx 
' r 2 " " x 
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Figure If-2. Geometric Parameters Used to Effect Transformation 
from Equation (II- 7) to Equation (II-8) 



H-7 



Thus, as programmed for the computer, the integral in question becomes 


2 r* 

•H-f W 


2 dx 
3 

xq 


L1-+XC+ 


i.iia-c 2 )] 2 E(k v 


... ;x2W:c2) Il-tec+^q ; x2)a : c2 


k (l^T 


.+4(1-x 2 )(1-c 2 ) 


k 4 (l-k 2 ) 


(n-8) 


The parameters q and k were defined earlier, and c = cos 0 2 - It should be ex- 
plicitly noted that the x in the denominator of the integral in Equation (11-8) would 
make the integral diverge, if the lower limit would go to zero. As will be explained 
in the next few paragraphs, there is good physical reason to put a finite lower 
limit at E. Unfortunately, however, the integral depends (by factors of 2 or 3) on 
just where the lower limit is placed. 


b. The Lower Limit Question 

Ligure 11-3 indicates the maximum that 6 ^ can be when considering 
the spherical character of the planetary surface. This largest value , ^im> satisfies 
the relation: 


lm a+h 


(II- 9) 
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Figure II-3. Geometry for Estimating E = Cos B 



where h is the altitude of the planetary probe, and a is the radius of the planet. 
From Equation ( 11 - 9 ), it should be clear than sin is near unity and, therefore, 
E = cos0i m is a very small quantity. From this, we can write: 


sin© 


lm 



1 2 , h 

— CZ HI--. 


(II-l 0) 


€ ^ 



This relation sets an estimate to the lowest value that the lower limit of the in- 
tegral in Equation (Pl-SJcan have. Unfortunately, for sufficiently small h, the 
use of this estimate for E yields physically unacceptable results; namely, that 
the average scattered power can be larger than the direct power. However, 
further thought suggests some way out of the difficulty. The lower limit , E arises 
from that part of the scattering surface which is very near .grazing. But that 
part of the surface would be very substantially shadowed by the roughness. One 
might, therefore, postulate rather reasonably that only part of the surface which 
receives rays at angles above the rms slope of the rough surface can contribute 
to the average scattered power. Conversely, this means that part of the surface 
which is near the grazing angle and receives rays from angles less than the rms 
slope is so heavily shadowed that very little energy is received (or scattered) there. 


With the aid of Figure 11-4 and some geometrical considerations, a reasonable 
estimate can be made for e using the above considerations, hvom Figure 11-4 
and the law of sines : 
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Figure II— 4. 
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Geometry for Estimating Shadowing Near Grazing Angles 
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or 



sin(s+p= sin0 lm 


a + h. a 

cos s X _2~ 

1+h/a ¥ lm 

j 7 

= v§ 2 * 2b/ft 


(n-ii) 


This particular formula suggests that there exists a lower limit to e even as h 0. 
From radar measurements of the lunar and Martian surface, it has been estimated 
that s lies between about 4° and 7°. Chdosing s = 0. i («5. 7°), Ewas calculated 
from Equation (Il-ii) , and that value of e was inserted into Equation (n-8). The 
results from numerical integration of Equation (n-8) are plotted in Figures H-5, - 
6, and -7, where, for purposes of estimating the Fresnel reflection coefficient of 
the Martian surface, it was assumed that E = 3 and a = 0. 


2, Interpretation of Results 

It must be admitted that it is very difficult to draw a firm conclusion 
from the results presented in Figures 11-5 through II- 7. hr fact, it appears that 
there exists a fundamental defect in the scattering formula, a defect which makes 
the whole computation for average scattered power quite suspect. It was decided 
to study this problem in greater detail. Reference is made to Equations {II— i ) and 
(n-3) and Figure II- i. In this report, P 0 dS represents the average power scattered 
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(b) HORIZONTAL POLARIZATION (COS 0 2 ! 


Figure n-5. The (Results of Integrating Equation 11-8 to Obtain 
the Average of Scattered Power for p = 0.1 




(a) VERTICAL POLARIZATION (COS e 2 ) 



(b) HORIZONTAL POLARIZATION (COS d 2 ) 


Figure 11 - 6 . The Results of Integrating Equation II— 8 to Obtain 
the Average Scattered Power for p = 1.0 
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(b) HORIZONTAL POLARIZATION (COS 0 2 ) 

Figure n-7„ The Exults of Integrating Equation H 8 to Obtain 
the Average Scattered Power for p = 10„ 0 
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by an element of surface, dS, divided by the power in the direct beam. The direct 
beam has a power flux density Pd ~ ^ > while the beam incident on the surface 
element dS is Therefore, the power scattered by a perfectly reflecting 

surface element, dS, relative to the power in the incident beam is 


p o r ds ■ ^ 

1 r 2 


(IS -12) 


2 [l+cos0 1 cos© ( ,-sin9 1 sin0„coscp] / 

F «W<» - ^ 


(cOS0,+CO [1 + p 1^] 


3/2 


and f 2 has already been defined in Equation (11-3). From conservation of power 
considerations, a requirement must be imposed on the function F; namely, 


dS 


jj F(0 1 ,0 2 ,cp) 


sin02d02 < hP = cos0^ dS 


(11-13) 


Equation, (n-13) is the mathematical statement that the total power scattered over 
all directions by a perfectly reflecting surface element [given by the left side of 
Equation (11-13)] must be equal to the incident power intercepted by the surface 
element dS [given by the right side of Equation (n-13)]. In addition to the con- 
servation of power requirement specified by Equation (n-3) , reciprocity con- 
siderations impose another important constraint; namely, ‘that the function, F, be 
an even function of <p and symmetric in 0* and 0 2 [i.e. . F( 6^, 6%, cp) = F (0 2 > -<p)]. 
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This constraint imposed by reciprocity means that if Equation (11-13) is not auto- 
matically satisfied, it is probably not possible (except fortuitously) to find a 
T i fudge factor" which would make the resulting scattering function satisfactory 
from both points of view- symmetry and conservation of power. 

The integral on the left side of Equation (11-13) was evaluated with the aid of a 
computer and is identified as Figure 11-8 presents a comparison of Nfy) 

with the required result, eos^. It can be seen that for p ^ 1, 0, N(0 1 ) does not 
look like cos at all. For p = 10.0, the effective scattering region tends to be 
reasonably concentrated around the specular point, and N(0^) begins to have some 
general resemblance to cos^. Out of desperation, it was decided to modify the 
previous F(0 1 , 8 2 ,<p) by multiplying F by a "fudge factor" cos Q L cos0 2 . This 
leaves the scattering function symmetric, but the requirement imposed by Equation 
(n-13) now becomes 


cos 0 2 F(0 1 ,0 2 ,cp) sin0 2 d0 2 dcp = 1 


(TI-14) 


This integral was also evaluated on a computer and is presented in Figure if-9as 

M(0 1 ). It can be seen that M(0 1 ) is not a constant as required. One must conclude, 
therefore, that the diffuse scattering function used in this analysis is not satisfactory 
for evaluating the average scattered power (or the fading rates) from an extended 
area. There is a fundamental defect in the scattering function unless p>10, and 
no reasonable amount of "fudging" can doctor it up. Perhaps this is because there 
is an inherent inconsistency in the whole theory when an exponential correlation 
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Figure II S. A Comparison of N (60 '.vith 
Cos Qi (From Equation 11-13) 
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function is assumed; namely, when the scattering function is derived based on a 
formulation for single scattering, the final result is a diffuse scattering function, 
which implies multiple scattering. 

It is of importance to note that the earlier study on scattering’ involving as it 
does highly directive scattering does satisfy both symmetry and conservation of 
power requirements as specified for Equation (11-13). 
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SECTION 111 


EXPERIMENTAL INVESTIGATION OF FAST FADING 

A . INTRODUCTION AND SUMMARY 

The experiment described in this task is an extension of the work done in 

an earlier study,' "Selected Studies of VHF/UHF Communications for Planetary 

\ 

(Mars, Venus) Relay Links." In that study, an experiment was conducted to 
evaluate the performance of a wideband noncoherent binary FSK system under 
the following conditions: 

© No fading in the presence of Gaussian noise; 

® Slow Rician fading in the presence of noise; and 
© Hardwire synchronization. 

The objective of this task is to extend the perforniance evaluation for conditions 
of fast Rician fading with and without doppler offset. The mathematical model 
of the channel is identical to that in the earlier study and is shown in Figure 
111-L For this evaluation, the Rayleigh generator was operated at fading rates 
greater than the bit rate (i.e., the correlation of the Rayleigh- distributed signal 
becomes very small over a bit interval) and the reflected- and direct- signal 
center frequencies have been offset to simulate differential doppler. 
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Figure 1114. Mathematical Representation of Ames Probe FSK Channel 











The following basic results were obtained: 

1. The effect of a reflected signal with a fading bandwidth smaller than 
the bit rate (slow fading) is to degrade the link performance relative 
to the no-multipath case. In terms of the mission, this means that if 
there are reflections from the planetary surface and the time-rate 
variation of the received reflected signal is slow relative to the bit 
time the link will be degraded. 

2. The effect of increasing the fading bandwidth (for a constant reflected 
power) is to improve the link performance. The performance for fast 
fading (large fading bandwidth relative to bit rate) is superior to the 
no-multipath case. In terms of the mission, this means that fast fading 
does not degrade the link. However, the fact that the fast-fading signal 
may actually improve the performance cannot be utilized since at this 
time, the existence of a given magnitude signal can only be postulated, 
but not guaranteed. 

3. The effect of an increasing frequency difference (differential doppler) 
between the reflected and direct signals is to improve link performance; 
the improvement is similar to that achieved by increasing the fading 
bandwidth. The impact of a differential doppler larger than the data 
rate is that the existence of multipath will not degrade the link. Again, 
since the existence of a large reflected signal can not be guaranteed, 

an improvement in n tk performance due to increasing doppler offset 


cannot be utilized. 



4. The limit of improvement in link performance, due to fast fading, 

increasing doppler offset, or any combination of the two, is to consider 
the received signal energy as the sum of the direct and reflected energy. 

B. EXPERIMENT IMPLEMENTATION 

The equipment used in this experiment is similar to that employed in the 
original investigation? The only differences are the use of a fast-fading gener- 
ator with provision for offsetting the reflected- signal center frequency from the 
direct- signal center frequency in place of the slow-fading generator. The ex- 
perimental system block diagram, shown in Figure III-2, is the same for both 
experiments. A block diagram of the fast-fading generator is shown in Figure 
in-3. 



Figure IH-2, Experimental System Block Diagram 

The fast-fading generator operates on the principle that the envelope of a narrow- 
band noise signal centered about a carrier frequency is Rayleigh distributed 2 . 

To implement this principle, a wideband noise source was operated into a narrow- 
band filter centered at 1 00 kHz. The physical implementation consisted of a 
General Radio noise source and a He wJcil -Packer d Model 310 Wave Analyzer. 
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Figure IH-3. Fas Fading Generator, Block: uiagram 
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WITH DOPPLER OFFSET 











The wave analyzer had selectable bandwidths of 200, 1000, and 3000 Hz. This 
provided a variable-fading-bandwidth capability, 

The signal at the output of the narrow-band filter can be represented by 

S R (t) = B* sinct^t+ B 2 cos aj R t, (IH-i) 

where atg is the center frequency of the bandpass filter (in radians per second). 

Bj 2 = B 2 2 = P R <m-2) 

where 

Bj 2 is the average value of the square of B l5 and 
P R is the average power of signal S R (t) . 

Then, S R (t) represents the reflected signal from the planetary surface with a 
power P R and a carrier frequency cn R . The envelope of S R (t) is Rayleigh 
distributed. 

(Ill- 3) 

(III-4) 


The drrect srgnal can be represented as 


Sj) (t) = A sin c<j£)t 

where the power of the direct signal is 
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Then, the sum of S R (t) + Sg(t) corresponds to the total signal received, or 

2 Sg(t) + Sg(t) = A sin togt + _B 2 cos togt + B t sin togt (HI-5) 
We can write 

W R = 0J D + a; B (m-6) 

where tog corresponds to the offset frequency between the direct and reflected 

signals. 

Now, 

A 

2 Sg(t) + Sg(t) = A sin togt + Bj sin (tog + tog) t + 

B 2 cos {tog + tog) t 

= [A + B* cos togt +B 2 sin togt] sin togt + 

[B 2 cos togt - B t sin togt] cos Wgt ( III— 7 ) 

The envelope of this expression is 

E (t) = 4 A 2 + B* 2 + B 2 2 + 2 ABi cos togt + 2 AB, sin togt (III- 8) 

For the case where there is no offset, tog equals zero; then 

E 0 (t) = VA 2 + B* 2 + B 2 2 + 2 AB X (IH-9) 

The envelope distribution tg Rician 3 where the reflected power is considered as 
noise power in Reference 3. 



The ratio of direct to reflected power (y 2 ) is from Equations III-2 and HI-4. 


2 _ Pp _ A 2 _ A 2 
r 2Bi 2 Bj 2 + B 2 2 


(m-io) 


The distribution is discussed by Schwartz 4 as a function of y 2 , where in his dis- 
cussion, y 2 would be the signal-to-noise power ratio. Where y 2 = 0, the 
distribution becomes Rayleigh and asy 2 >> 1, the distribution approximates a 
Gaussian distribution about the amplitude A with a variance Bj 2 . 


The output of the fast-fading generator for several values of fading bandwidths 
and •'} is shown in Figure III— 4. These photographs were taken from an oscillo- 
scope connected to the output of the bandpass filter shown in Figure III-3a. The 
wave analyzer (bandpass filter) was operated in the automatic-frequency-control 
mode which locks the center of the bandpass to the frequency at the input, thus 
ensuring that the offset between the direct and reflected signals is zero. The 
signal at this point is not frequency shift keyed (FSK) . The direct component 
power was held constant in all photographs except fory 2 = 0, where the re- 
flected power was made equal to the direct power in the other photos (direct 
power is zero fory 2 =.0) . The vertical and horizontal scales are constant 
throughout the sequence, 


Observing any horizontal set of Figure IH-4, it is seen that the amplitude dis- 
tribution is about the same for all three bandwidths but that the time distribution 
varies with bandwidth. This is expected, since Equation III— 9 shows that the 
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DIR. PWR/REF. PWR FADING BANDWIDTH (0J B =O) 
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HORIZONTAL SWEEP RATE = 5 MILLISEC/ctg 
VERTICAL SENSITIVITY = 0.2 V/cm 










amplitude distribution is only a function of the powers of the direct and reflected 
signals and not of their bandwidths. An important aspect to note is that for the 
experiments conducted here, the bit time was 10 milliseconds and that as the 
bandwidth is increased, a larger number of amplitude variations occur within a 
bit period. This means that if the envelope amplitude is averaged over the bit 
period, more independent samples of the waveform are obtained for the larger 
bandwidths and the measured average value tends to approach the true average 
value with a higher probability. In mathematical terms, 

T 

^ J E (t) dt “ E asB>>-^ (111-ii) 

J o 

where E is the average value of E (t), B'is the fading bandwidth, and T is the 
bit period. 

Observing the vertical sets of Figure IU~ 4, as y 2 increases, the magnitude of 
the amplitude variations tend to decrease. Wheny 2 >> i, the distribution tends 
to become Gaussian about the direct- signal amplitude, A, and the variance 
approaches the square root of the reflected power, /b^ 2 . 

The effect of frequency offset between the direct and reflected signals is shown 
in Figure III— 5 for y 2 = i, and B = 200 Hz. The basic equation describing the 
envelope for this case is Equation III- 8. 

E (t) = 7a~+ Bj 2 *- B 2 2 + 2A [Bj cos ce-gt-.- B^, sin Wgt] 
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Figure 111-5. Fading Simulator Output Waveforms (with Doppler Offset) 
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VERTICAL SENSITIVITY = 0.20 V/cm 



When a>n >> B (the fading bandwidth), the term Bj cos OJjjt + B 2 sin q>q t becomes 
a narrow-band Gaussian distribution centered about ojg (B t and B 2 are each 
independent Gaussian variables) and can be written as 

cos Wjjt + B 2 sin c^t =* Q sin (otgt + Q ) , (In-12) 

where Q = ^ B,’ -*■ (Rayleigh distributed), and/) is a random variable with 

all values equally likely for 0< B< 2 tt, Thus, Equation D3-8 can be written as 

E = V A 2 + Q 2 + 2 AQ sin (Wgt + 9 ) - T (111-13) 

The 2AQ sin + 9) component can readily be seen in Figure III-5 for ojg = 

5Q0 Hz and OJg = 1000 Hz. For OJ-g = 3000 Hz, the horizontal time scale does 
not allow resolution of the offset frequency. If only the envelope of the 2AQ sin 
(Wgt + 6) component is considered (as in fact occurs in the 0)g = 3000 Hz photo- 
graph of Figure 111-5) then the overall envelope is 

E = VA 2 + Q 2 + 2AQ or E = A + Q (111-14) 

which is a constant (A) plus a Rayleigh variable (Q) . By observing the CUg = 3000 
Hz recording in Figure 111-5, the appearance is as predicted. A Rayleigh wave- 
form is shown for comparison in Figure HI-4 (y 2 = 0) without the constant A. 

A detailed block diagram of the equipment implementation employed in the tests 
is shown in Figure III— ^ , ™i'th. an equipment chart. Thf 1 FSK receiver is the same 
as employed in the original study' and is described in detail in the final report. 



Figure III-6. Equipment Implementation Block Diagram 
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C. EXPERIMENTAL RESULTS 


1. General 

The results of the experiment determine the effects of fast fading (with 
and without doppler offset of the reflected signal) on a wideband noncoherent FSK 
system. The mathematical analysis in the previous paragraph is now continued to 
determine the effects of slow fading and fast fading on the FSK system. The 
analysis is carried out sufficiently to qualitatively describe the effects of fading on 
the FSK system as a function of the pertinent fading parameters, i. e, , fading band- 
width, ratio of direct signal power to reflected power (y 2 ), and doppler offset. The 
mathematical analysis is presented in order to provide a background to evaluate 
the experiment, 

2, Mathematical Analysis 

An assumption implicit in all the analysis and experimental data is that the 
time delay between a bit arriving at the receiver along the direct path and one 
along the reflected path is negligible compared to the bit period. The low bit 
rates and the flight geometry for the missions considered justify this assumption. 

The envelope of the incoming signal is described in Equation HI— 8 as 

E(t) = / A 2 + Bi 2 + B 2 2 + 2A B x cos U> B t + 2AB 2 sin 

For the case where WT is large, it was shown in Reference I (Section IH (E)) 
that the signal output of the envelope detector is properties 1 to the square of 
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its input and that the noise power output is a constant. Therefore, the- signal 
at the output of the detector can be written as 


E*(t) - K(A 2 + Bi 2 H-Bg 2 + 2AB* cos Wgt + 2AB 2 sin Wgt) 
where K is a detector constant. 


(In-15) 


The integrate and dump circuit (matched filter) performs the operation of averag- 
ing the signal over a bit period so that 


E* = — j / a 2 a. r>2 j. -d 2 

T I 

Jo 


- J (A 2 + B t 2 + B 2 2 + 2AB X cos Wgt + 2AB 2 sin Wgtjdt (IH-16) 

o 


The noise output of the same circuit is 


N* 


■if 


Ndt 


(In-17) 


where N is a Gaussian distributed noise variable. The output signal-to-noise 
ratio is then 


mmm /•T 2 2 2 

M- = — j (A +Bi +B 2 +2AB 1 cos W B t+2AB 2 sin W t) dt (UI-18) 
N 1 T / 


L 


Ndt 


This is a general expression of performance as a function of fading bandwidth 
including doppler offset. 


In order to evaluate the effects of different fading conditions, a parameter, p, 
is defined, which is the ratio of the signal-to-noise ratio in :he link with 
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multipath present to the signal-to-noise ratio without multipath (all other con- 
ditions remaining constant). From Equation ni-18, the signal-to-noise ratio 
with no multipath (B* = B 2 = 0) is 


S 

N 


1 rT. 

T *f) 

i r T 

T * / o 


A 2 dt 
Ndt 


1 

T 



(m-i9) 


The ratio of S/N’s with and without multipath from Equations UJ-18 and 111-19 is 
then 


B l 2 + V B 1 ®2 

p = _ (1+ — - — ' + 2 ~t cos Ct) t + 2 ~ sin Wt) dt (HI-20) 

i . w A d A 

A 


This is the generalized expression for p and the intent is to analyze it for the 
various cases of slow fading and fast fading with and without offset in order to 
provide a theoretical background for the experimental data obtained. 


a. Case I: Slow Fading, No Doppler Offset 

This is an experiment that was performed during the earlier study 1 
and is reintroduced here in summary form for completeness and to emphasize 
the contrast between slow fading and fast fading. 


Bi and B 2 are random variables with Gaussian statistics. For slow fading, B* 
and B 2 can be considered constant during the bit period since BT «1. 



Therefore, 


B„ 


B, 


7-f 


dt 


*T 2 
B 2 dl 

7 d ‘ 


(m-2i) 




0 (nodoppler offset) 


For this case, p (from Equation III-20) reduces to 

2 2 

B 1 +B 2 2B 1 
A = 1+ ^? + 1T 


(in-22) 


and 


2B 1 

A 


|> / ~ ^2 \ , /; will be less than one and the performance of the link 

A / 


"A" + 


b. Case II; Fast Fading, No Doppler Offset 

In the case of very fast fading, B]_ and Bg change rapidly over a bit period. 

Since Bi and B2 are Gaussian-distributed random variables, and their average value 
is zero, then 



CALCULATED ERROR RATE 
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f 1 f* 

A / B 1 dt — O A / 

*X> Jq 
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B 2 2 dt 


E„2 (HI-24) 

* A 


0) = 0 (no doppler offset) 

D 


For this case, p (from Equation m-20) reduces to 


= 1 + 


- 2- 2 
B 1 +B 2 


(m-25) 


n 


The ratio of direct signal power to reflected power (y2 ) is equal to — g 

®1 +S 2 

V 

therefore, p can be written as 


n 


= l + 


( ni-26 ) 


For fast fading, one would expect the link performance to be improved by the ratio 
Pjj above. The improvement approaches PjjWhenBT approaches infinity, where 
B is fading bandwidth and T is bit period. 

The experimental data for a fading bandwidth of 3,000 Hz are shown in Figure 111-8. 
Since the data rate for these experiments is 100 bits per second, a 3000-Hz fading 
bandwidth corresponds to BT of 30. Aty2 = there should be an improvement over 
y2 =ro of 3 db, andaty2=4, an improvement of ldb. For .amean error rate (Pe) 
of 10“3 ? these predictions are realized, At error rates of 10"^ and 10~5 j the cor- 
respondence between prediction and experiment is not as good indicating that the 
fading bandwidth, B, needs to be larger to ensure that B^ dt is small 

enough. The design center for this mission is a mean error rate of 10 -3 . 



BIT ERROR RATE (P e ) 



Figure 111-8. Hit Error Rate as a Function of Signal -to-Noise 
Ratio and y z for a Fading Bandwidui of 3 kHz 
(No Doppler Offset) 




The experimental data for fading bandwidths of 200 Hz and 1000 Hz, which cor- 
respond to BT’s of 2 and 10, respectively, are shown in Figures 111-9 and III-1Q, 
These curves show the degradation of link performance as B becomes small. 
Even though BT is small for the 200-Kz fading bandwidth, the performance is 
still far superior to the very slow fading case cf Figure 111-7. The effect of BT 
(fading rate) on the error curves is indicated in Figure III— 11 , where a family 
cf curves for various BT ! s are plotted for y 2 =? 1. 

A set of data was taken for a bit rate of 1000 bits-per second and a fading band- 
width of 1000 Hz„ This corresponds to a WT product of 10 and BT of 1, These 
curves are shown in Figure III-12. All previous data presented were taken at a 
WT product cf 100. The bit error rate performance at WT = 10 is consistent 
with that taken at WT = 1000 and the same general trends are exhibited, as 
would be expected from the analysis. 

The application cf the data on fading to the mission requirements is perhaps 
best presented in the form cf the penalty imposed by the presence of a reflected 
signal (multipath) as a function cf BT. The mission design center is an error 
rate cf 10 -3 o Figure III-13 presents the margin required over an ideal link with 
no multipath to achieve an error rate better than 10“ 3 . From the previous 
multipath nnalysis, y 2 > 4 can be achieved with an omnidirectional antenna 
aboard the probe by constraining the look angle from the probe to the spacecraft 
bus from the zenith to 62 degrees. Additionally, the slowest fading bandwidth, 
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SIGNAL-TO-NOISE RATIO (DB) 


Figure 111 - 10 . Bit Error' Rate as a Function of Signat-to-Noise 
Ratio and y 2 for a Fading Bandwidth of I kHz 
(No Doppler Offset) 
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SIGNAL -TO- NOISE RATIO (DB) 

Figure ffl-li. Bit Error Rate as a Function of Direct Signal-to-Noise 
Ratio with Ratio of Fading Bandwidth to Data. Rate (BT) 
as a Parameter 
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Figure III-12, Bit Error Rate as a Function of Signal-to-Noise Ratio 
and y 2 for a Fading Bandwidth of 1 kHz and Bit Rate of 
i kHz (No Doppler Offset) 
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-13. Fading Margin as a Function of 




which occurs at terminal velocity, is approximately 200 Hz and corresponds to 
BT = 2 for the 100 bit-per- second data rate. For these considerations, a fading 
margin of 2 db would be required. If BT is greater than 4 aty 2 = 4, there would 
be no penalty introduced by planetary surface reflections. It should be pointed 
out here that Figure IU-13 applies for a WT = 100 and if this parameter is 
changed, a new set of curves would be required to determine the fading margin. 

It is for this reason that point A on Figure 111-13, which was taken from the 
WT = 10 data (Figure 111-12) fory 2 = 1, BT = 1, does not fit on the curve; 

An interesting phenomena illustrated in Figure 111-13 is that for large BT pro- 
ducts, the reflected signal aids the link performance. This is due to the fact 
that the FSK detector is basically an energy detector and if the reflected signal 
has zero correlation with the direct signal (Equation 111-21) over the bit period, 
the detector will see the sum of the engeries {Equation 111-25). The effect of 
rapidly varying scattering (fast fading) by a rough planetary surface is to convert 
the incident wave into noise-like reflected signal whose correlation with the 
direct signal decays rapidly with time. 

c. Case 111: Doppler Offset Larger Than Fading Bandwidth 
In this case, Wg »Band, therefore, 

B* cos cu-gt + B 2 sin cu-gt =Q sin (u>gt + 0) (In-27) 

where Q = VBj 2 + B 2 2 

and 0 can take on any value between 0 and 2 7T with equal iik'iinood. 
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Therefore, Equation 111-20 becomes 


= -L f?(l + T5 + ^ 2 ) dt + _2_ X T Q sin (w B t + 0) dt (IH-28) 
'In. T 0 *A2 A2 AT 0 

When many cycles of are included in a bit period (T), then 

~r f Q T B sin (Wgt + 0 ) dt — 0 as cOg — (In-29) 

Thus, for large Wgt products (Wgt^$>l): 

P TTT is always greater than unity, but is time varying as a function of Bi + B 2 . 


The results of an experiment with fading bandwidth of 200 Hz and variable dop- 
pler offset fory 2 = lare shown in Figure HI-14. It can be seen that with no 
doppler offset and in the region cf interest (P e = 10" 3 ), the performance is 
degraded with respect to the no-multipath case. However, as doppler offset 
is introduced, the performance improves rapidly to the case where the perform- 
ance with multipath is superior to that with no multipath Figure IH-15 is a 
similar set of curves for a fading bandwidth of 1000 Hz. In this case, the fading 
bandwidth was -such that even without doppler offset, the performance was 
superior to the no-miltipath case (y 2 = °°). The performance here improved 
with doppler offset to the limiting value predicted by 
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SIGNAL -TO-NOISE RATIO (DB) 


Figure HI-14,, Bit Error Rate as a Function of Signal-to-Noise Ratio 
with Variable Doppler Offset and a Fading Bandwidth 
of 200 Hz 











Since y 2 = l,the expected curve would be shifted 3 db to the left of the no- 
muitipath case. In this case, due to the high fading bandwidth, 

/ 0 T Bj 2 - Bl 2 , and ~ f Q T B* 2 - B 2 2 (IH-32) 

where bJ* = B. is the power in the reflected signal. 

Figure UI-16 shows a set of data for a 200-Hz fading bandwidth with 
3-kHz doppler offset and variable y 2 . The doppler offset is sufficiently great 
so that the performance is always better than for no multipath. However, the 
curves are not displaced by a constant 1 + ”ras would be the case for a higher 
fading bandwidth. The 1000-Hz offset curve for 1000-Hz fading bandwidth of 
Figure 111-15 does show a nearly constant offset. Even for a 3000-Hz fading 
bandwidth with no doppler offset (Figure 111-8), the curves are not as parallel 
as for the 1000-Hz offset, 1000-Hz fading bandwidth curves of Figure 111-15. 
Equation 111-8 shows that for sufficiently large fading bandwidth, the curves will 
be parallel without doppler offset. However, this ideal can be approached more 
rapidly by introducing some doppler offset. Therefore, the most desirable 
situation is to have a combination of fading bandwidth and doppler offset. For 
example, the data shows (Figure IU-14) for a fading bandwidth of 200 Hz, and a 
doppler offset of 200 Hz (the estimated terminal impact values), there would be 
performance improvement at the design center of P e = 10” 3 . This applies to a 
bit rate of 100 bits per second; however, the conclusion can be generalized by 
referring to frequency time products. 
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Figure HI-16,, Bit Error Rate as a Function of Signal-to-Noise Ratio 
with. Variable y 2 and a Fading Bandwidth of 200 Hz and 
a Doppler Offset of 3 kHz. ....... .. . 





Therefore, in general terms, if 


BT (fading bandwidth X bit period) >2 
and OU 3 T (offset frequency X bit period) >2. 

there will be no degradation of performance due to multipath and in fact, there 
will be some improvement depending ony 2 . The exact shape of the performance 
curves is determined also by the WT parameter and therefore, the numerical 
values presented apply only for WT = 100. However, the general conclusion still 
applies that the most effective mechanism to reduce the adverse effect of multi- 
path is to have moderate fading bandwidth with some doppler offset. 

Thus for the phenomena of link performance improvement due to multipath has 
been discussed. The purpose of the study was to determine the effect of a 
possible multipath environment to ensure the success of the link. At our present 
state of knowledge, we cannot guarantee a very specific set of surface charac- 
teristics in order to define the multipath environment precisely. The best that 
can be done with the limited knowledge available is to postulate a set of typical- 
case conditions to ensure adequacy of the link in the event that surface reflections 
are present but in no way count on them to be present to aid the link performance. 
Thus, multipath becomes a possible degradation factor but should not be used to 
aid link performance. 
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SECTIONIV 


EFFECT OF A BIT SYNCHRONIZER ON 
LINK PERFORMANCE (TASK 4) 

A. INTRODUCTION AND SUMMARY 

The previous study' experimentally evaluated the effects of slow fading 
(BT << 1 3>n the performance of the link with a bit synchronizer. The purpose 
of this experiment is to determine the effects of fast fading (BT — 2 ), with and 
without doppler offset, on the performance of the link with a bit synchronizer. 

In the tests, the FSK receiver was operated with bit synchronization derived 
from the data signal, as contrasted to Task 3 (Section HI) where the bit 
synchronization was hardwire-connected from the data source. 

The results of the earlier study showed a degradation of the order of 2 db at 
error rates of 10 -3 due to use of the bit synchronizer in the presence of slow 
fading (y 2 = 4 and y 2 = 10) . The degradation was even greater for larger error 
rates. The results of the fast-fading experiment showed no degradation due to 
the fading. A uniform degradation of the order of 0.5 db was introduced by the 
bit synchronizer; however, this degradation was noted whether or not there was 
fading present. 

B. EXPERIMENT IMPLEMENTATION 

The equipment used Cor this experiment is identical to that employeh in the 
slow-fading tests and is described in Reference 1, Section 5.2. The bit syn- 

IV- 1 

'i • - ■ 



chronizer employed was a Telemetries Model 6203A. The tests were conducted 
on "A", M B" comparison basis where at each signal-to-noise ratio, the bit 
error rate was first measured with hardwire sync and then measured with the 
bit synchronizer in place of the hardwire sync. 


C. EXPERIMENTAL RESULTS 

The effect of the bit synchronizer was first measured in the absence of 
fading (y 2 = °°) . The data are presented in Figure IV-I. The measured data 
points are shown plotted along with the theoretical (perfect sync) curve. The 
data show a small degradation, on the order of 0.5 db, for error rates between 
iCf 2 and 10 -3 . These results are consistent with those obtained in the previous 
study. 


The fast-fading tests were performed under conditions corresponding to three 
points in the sample mission, i.e.: 


Fading bandwidth = 3,000 Hz 
Doppler offset = 3,000 Hz 

Fading bandwidth = i, 000 Hz 
Doppler offset = i, 000 Hz 

Fading bandwidth = 200 Hz 
Doppler offset = 500' Hz 


(Just emerging from blackout) 

(Intermediate between blackout and 
terminal impact) 

(Just prior to impact) 
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The results of these tests are shown in Figures IV- 2, IV-3, and IV-4. The data 
show consistent displacement of the bit synchronizer curves of about 0.5 db; 
however, within the limits of experimental accuracy this degradation is no more 
than that observed without multipath. 

I*" 

The effect of using a bit synchronizer in the presence of slow fading is shown in 
Figure IV-5 to illustrate the pronounced impact of slow fading on link perform- 
ance. In the case of very slow fading, fades occur of sufficient time duration 
to cause loss of lock in the bit synchronizer. The reacquisition process then 
introduces a loss of data beyond that caused by the fading itself, thus in effect, 
increasing the error rate. For the case investigated here, 200-Hz fading band- 
width and greater, the probability of a fade being of sufficient duration to cause 
loss of lock is extremely small and’therefore, the results obtained are not un- 
expected. These results are further substantiated by the findings of Task 6 
(see Section V) which show that even for 200-Hz fading bandwidth without 
doppler offset, there is no discernable correlation between bit errors, or in 
other words, given that an error occurs, the probability that succeeding bits 
will be in error is not affected. This means there is no clustering of errors and 
the probability of a synchronizer loop dropout is small. 

The conclusion then, is, that the use of a bit synchronizer in the presence of 
fast fading does not cause system error rate performance degradation beyond 
that caused by the bit synchronizer in the absence of fadirg. 


IV-4 



BIT ERROR RATE (P e ) 



Figure IV-2. Effect V a pit Synchronizer on Bit Error Rate for a Fading 

Bandwidth of 3 kHz and a Doppler Offset of 3 kHz, y 2 Variable 
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Figure TV- 3 . 


Effect of a Bit Synchronizer on Bit Error Rate for a Fading 
Bandwidth of 1kHz and a Doppler Offset of i kHz, y 2 Variable 
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Figure IV-4. Effect of a Bit Synchronizer on Bit Error Kate for a Fading 

Bandwidth of 200 kHz and a Doppler Offset of 500 Hz, y 2 Variable 
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Figure N-5 . Effec" of a Bit Synchronizer for Very Slow Fading (BT « 1) 
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SECTION y 


LINK PERFORMANCE CRTTERLA (TASK 6) 

A. INTRODUCTION AND SUMMARY 

The performance of digital transmission systems is mos-t-often specified in 
terms of bit error rate. The assumption implicit in such a performance criteria 
is that the errors are random and independent. This means that the probability 
of an error for any given data bit is the same as for any other bit and that past 
errors and future errors do not influence the probability of error. In communi- 
cations systems that operate under conditions where white Gaussian noise is 
added to a fixed signal, the occurrence of bit errors is random and independent. 
In a multipath environment, the signal amplitude is not fixed, but varies in some 
random manner before the addition of the Gaussian receiver noise. Where 
multipath causes slow fading (fadingperiod is greater than a bit period), deep 
fades occur which result in a high error rate over the duration of the fade. 

Since the fading period is greater than the bit period, the probability of error 
is not independent from bit to bit. System performance in such cases is often 
specified in the percentage of time the system operates with a low error rate. 

An example of this is the use of percentiles for expressing the performance of 
troposcatter systems. 

0 

The experiment performed here is to determine whether the FSK system pro- 
duces random, independ'"* 1 ! errors under the expected r ‘ouuitions of fast fading. 
The system was tested over a, range of fading bandwidths and doppler offsets for 
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small sample lengths (700 bits) and the error distribution compared to that 
predicted by theory, assuming random independent errors. The results showed 
good correspondence between prediction and measurement indicating that the 
errors are truly independent and random. The conclusion is that for the range 
of parameters investigated ( BT > 2 and? 2 >1 ), the average bif error rate is a 
valid measure of performance. 

B. EXPERIMENT IMPLEMENTATION 

The equipment used for this experiment was the same as that employed for 
implementing the experiment of Task 4. The error counter was programmed to 
read the number of errors in a 700- or 1400-bit sample rather than a 10,000- or 
100,000-bit sample. In Task 4 (SectionlV) , the objective was to obtain the 
average error rate; therefore, very long sample lengths were necessary to 
ensure a sufficient number of errors in the sample to provide a high confidence 
that the average had been approached. The purpose of this task is to observe 
the errors in many small samples to determine the distribution of the errors 
among the samples. 


C. EXPERIMENTAL RESULTS 

The probability of K errors occurring in a sample of 7] bits for independent, 
randomly occurring errors is given by: 


ic = C K ,?(P e )K <V 


(V- K) 


(V-i) 


V-2 
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where 


c k t? 

p e 

Qe 


is the probability of K errors; 

is the combination of 77 things taken K at a time; 

is the probability of a bit error; and 

is the probability of no error = (i-Pe) . 


Equation V-i is the binomial distribution which can be written as 


P K = 


77! (P e ) K (1-P e ) (?7 ~ K) 
K ! (77 -K) ! 


(V-2) 


The average or expected value of K is 77 P e and the variance is 77 P e (i-P e ) . 
When 77 is large compared to K ( 77 »K) and P e is small compared to unity, 
the binomial distribution is approximated by the Poisson distribution, 


p = (77 P e ) K e" 7 ? P e 
K K! 


(V-3) 


In the present experiment, the Poisson distribution was used because the 
approximation is very good since P e is on the order of iO -3 and K is less than 
10 compared to an 77 of 700. For the case P e <<l, the expected value and 
variance of the binomial distribution is the same as that for the Poisson 
distribution. 


From the experimental data, the total number of bits and total number of errors 
were calculated and tha probability of error for each run was computed by 
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n 


p e - 


m 


( V-4) 


where 

nQ is the total number of errors; 

N is the number of samples (approximately 100 samples per run) ; and 
77 is the number of bits per sample. 


This then gave the average probability of error for the overall run, The 
measured P e was used in the Poisson equation to predict the error distribution. 
The number of samples containing K errors was then plotted as a function of 
K to give an error- density distribution. The predicted values were then com- 
pared to the measured values in graphical format. 


The tests were performed over a range of conditions similar to that which would 

be experienced during a typical mission profile. The intention was to perform 
the tests at a signal-to-noise ratio corresponding to a bit error rate of io~ 3 . 

The mean bit error rate is a calculated value from a run of data. In order to 
approach the desired bit error rate, a lengthy reiterative procedure is required 
(approximately 30 minutes per iteration or run). The bit error rate is quite 
sensitive to small dhanges in input signal-to-noise ratio which is variable in 
one-db steps. In most cases, a one-db change in signal-to-noise ratio pro- 
duces an order of magnitude change in bit error rate in the vicinity of interest. 
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(The experimental curves of Task 3, for fast fading, demonstrate this character- 
istic.) Therefore, although the objective was to perform the tests at P e = i0” 3 » 
the actual values of P@ varied from 5 x i0~ 4 to 2 x iO” 3 . This includes the effect 
of equipment stability, which is minor, during the runs. 

**■ 

The first two runs were performed with no fading; one for a sample length of 
700 bits and one for a sample length of 1400 bits. The results are presented in 
Figures V-l and V-2. The measured distribution of errors is in close agreement 
with the predicted distributions for independent, randomly occurring errors. 

This is expected for the non-fading case and indicated that the experimental set- 
up is operating satisfactorily. 


The next series of runs was performed under fast-fading conditions simulating 
three phases of the probe descending through the planetary atmosphere as follows: 


Phase I ~ Just After Emerging from Blackout . 
B = 3000 Hz y 2 = 10 Figure V-3 

COg = 3000 Hz y 2 = 4 Figure V-4 


Phase II ~ Midway in the Mission . 

B = 100 Hz y 2 10 Figure V-5 

COrj = 1000 Hz y 2 = 4 Figure V-6 


Phase in 
B = 200 Hz 
= 500 Hz 


~ Just Before Terminal Impact . 

y 2 = 10 Fig'-.xe V-7 

y 2 - 4 Figure V-8 
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Bit Sample 
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Figure V-3 Er*- 0 * Distribution for B = 3000 Hz, cu B = 3000 Hz with y 
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Figure V-5. Error Distribution for B = 1000 Hz, co B - 1000 Hz with y 
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Figure V-6 Distribution for B = 1000 Hz, w = 1000 Hz with y 
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Figure V- 7. Error Distribution for B = 200 Hz, o> = 500 Hz with y 
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Figure V-8. Error Distribution |r B = 200 Hz, co - 500 Hz with y 





In each of these tests, the measured results compared favorably with the pre- 
dicted values indicating that the multipath (under the stated conditions) did not 
cause a predisppsition towards multiple errors. These results are not unex- 
pected since as noted in Task 3 under the same conditions (Figures 111-14, 

IE- 15, and III-i6), multipath fading did not cause degradation of the link and 
actually provided improvement over the no-multipath case. 

An additional series of tests was performed with fading bandwidths of 200 Hz 
with no doppler offset for y 2 = 10, y 2 = 4, and y 2 = 1, Under these conditions 
in Task 3, there was pronounced degradation of the link due to the presence of 
multipath fading. If an effect on the error distribution is to occur, it would be 

more likely to occur under these conditions, since BT = 2, o> =0 and y 2 = 1 are 

B 

the expected worst cases. The results of these tests are shown in Figures V-9, 
V-iO, and V-il. Even under these conditions, there is no experimental evidence 
of a higher incidence of multiple errors. This indicates that there is very little 
correlation of the fading effects from one bit period to the next. These results 
were surprising in view of the performance impact. The waveform photographs 
of Figure 111-4* however, provided additional qualitative evidence to support 
the results. Observing the waveforms for the 200-Hz fading bandwidth over a 
bit period (10 milliseconds), there is no obvious correlation with the waveform 
in any other bit period. 
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v-15 


Figure V-9. Error Distribution for B = 200 Hz - 0 with y 2 - 10 
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Figure 11. Error Distribution for B = 200 Hz, Oio = 0 with 





D. CONCLUSIONS 


In the case of very slow fading, it is known that deep fades occur which 
cause a large increase in the probability of error during that period of time. 
Therefore, in sampling the data, a higher number of errors would be expected 
during a fade, while conversely, a very small number of errors would be ex- 
petted during periods of reinforcement (which are equally likely as fades). 
The error distribution would then show a disproportionate tendancy toward ex- 
tremes of small and large numbers of errors relative to that predicted by a 
mean error rate. The results of this experiment have shown that for BT pro- 
ducts greater than two, this effect has not occurred. The errors are random 
and independent, with no tendency to cluster, and the average bit error rate is 
a valid measure of communications performance. 



